Dynamics of entanglement for coherent excitonic states in a system of two coupled 

quantum dots and cavity QED 
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The dynamics of the entanglement for coherent excitonic states in the system of two coupled large 
semiconductor quantum dots (R/üb í§> 1) mediated by a single-mode cavity field is investigated. 
Maximally entangled coherent excitonic states can be generated by cavity field initially prepared 
in odd coherent state. The entanglement of the excitonic coherent states between two dots reaches 
maximum when no photon is detected in the cavity. The effects of the zero-temperature environment 
on the entanglement of excitonic coherent state are also studied using the concurrence for two 
subsystems of the excitons. 
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I. INTRODUCTION 

Demonstration of the existence of fast quantum algo- 
rithms [1] excited scientists to pour tremendous enthu- 
siasm into the quantum computation and information 
theory. In this relatively new area, scientists not only 
make abstract theory, but also make efforts to find some 
physical systems to test and realize their ideas. A rep- 
resentative list includes the trapped ions [2], the liquid- 
state nuclear magnètic resonance [3] and the solid-state 
material [4]. 

Nanotechnology opens technological possibilities to 
fabricate mesoscopic devices. Semiconductor nanos- 
tructures, especially quantum dot structures, are very 
promising for the realization of quantum computation 
and the quantum information processing. The quantum 
gate realization using quantum dots has been proposed 
[5]. Sanders et al. have discussed the scalable solid- 
state quantum computer based on elèctric dipole tran- 
sitions within coupled single-electron quantum dots [6]. 
Imamoglu et al. propose a scheme which realizes con- 
trollable interactions between two distant quantum dot 
spins by combining the cavity quantum electrodynamics 
(QED) with electrònic spin degrees in quantum dots [7]. 
Entanglement of the exciton states in a single quantum 
dot or in a quantum dot molecule has been also demon- 
strated experimentally [8,9]. The reference [10] investi- 
gate the entanglement of excitonic states in the system of 
the optically driven coupled quantum dots theoretically 
and propose a method to prepare maximally entangled 
Bell and Greenberger-Horne-Zeilinger states. Berman et 
al. has investigated the dynamics of entangled states and 
a quantum control-not gate for an ensemble of four-spin 
molècules and shown that entangled states can be gener- 
ated using a resonant interaction between the spin system 
and an electromagnètic 7r pulse [11]. 

The main focus of the references cited above is the en- 
tanglement of the excitonic orthogonal states. In fact, 
not only the entangled orthogonal states but also entan- 
gled nonorthogonal states play an important role in quan- 



tum cryptography [12] and quantum teleportation [13]. 
One kind of entangled nonorthogonal states is the en- 
tangled coherent states [14,15] which can be generated 
by a nonlinear Mach-Zender interferometer. However up 
to now, no scheme for generating entangled excitonic co- 
herent states has been proposed. This paper will focus 
on this tòpic, especially, on the dynamics of the entan- 
gled excitonic coherent states in the system of two large 
semiconductor quantum dots coupled by a single-mode 
field. The conditions for the preparation of maximally 
entangled excitonic states are derived assuming that the 
cffective radius R of each quantum dot is much larger 
than the Bohr-radius ag of the bulk-exciton. The inter- 
actions between two quantum dots is mediated only by 
the single-mode cavity field. 

This paper is organized as follows: In Sec. II, we will 
describe the Hamiltonian of the system of interest, and 
obtain the solutions of the operators which present dy- 
namical quantity. Then we will give the wave function of 
the system with certain initial states. In Sec. III, we will 
discuss the entanglement of the excitonic coherent states 
in the system of the two coupled quantum dots by virtue 
of the concurrence [16,17]. The effects of environment on 
concurrence will be discussed in Sec. IV, and finally the 
conclusion of this study will be presented in Sec. V. 



II. HAMILTONIAN AND SOLUTIONS 

The model that is analyzed in this study consists of 
two quantum dots which are placed into a single-mode 
cavity. The quantum dots have large sizes satisfying the 
condition R^> as- We also assume that there are a few 
electrons excited from valence-band to conduction-band. 
Then the excitation density of the Coulomb-correlated 
electron-hole pairs, excitons, in the ground state for each 
quantum dot is low. This, in turn, implies that the av- 
erage number of excitons is no more than one for an 
cffective area of the excitonic Bohr radius. Therefore 
exciton operators can be approximated with boson oper- 
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ators, and all nonlinear terms including exciton-exciton 
interactions and the phase space filling effect can be ne- 
glcctcd. The ground energy of the excitons in each dots is 
assumed to be the same. There is a resonant interaction 
between the single-mode cavity field and the excitons. 
The distance between two quantum dots is assumed to 
be much larger than the optical wavelength A of the cav- 
ity field, so the interaction between quantum dots can be 
ncglcctcd. Under the rotating wave approximation, the 
Hamiltonian for this system can be written as [18] 

2 2 

H = huià^à + hui^2 HJ>m + 9m{blnà + a f è m ) 

m— 1 m— 1 

(1) 

where (à) are the operators of the cavity field with fre- 
quency uj, &J„(6 m ) denote the exciton operators with the 
same frequency u> of the cavity field, and m — 1(2) rep- 
resent the first dot or the second dot. The coupling con- 
stants between quantum dot one (two) and cavity field 
are represented by g m with m = 1,2. Without loss of 
generality, we can take the two coupling constants be- 
tween the cavity field and the quantum dots as different. 
The Hcisenberg equations of motion for the operators of 
the cavity field a(a t ) and the excitons b m (w m ) can be 
easily obtained as 



— = -iaia - igxbi - ig 2 b 2 
ot 

dh . f . . 
— = -iwbi - igia 
ot 



dt 



-iub 2 — ig 2 à 



(2a) 
(2b) 
(2c) 



It is very easy to obtain the solutions of the above oper- 
ator equations as 



o(í) - cos( ff í)a(0)e~ Mt - i^- sin(.gí)6 1 (0) e -^ t 



The two modes of excitons are in the vacuum states 
|0)i|0)2 with label 1(2) denoting exciton mode one (two). 
The whole initial state for the excitons and the cavity 
field can be expressed as 



-i^l S in( 5 í)è 2 (0)e- ia,í 
6i(í) = -i— sm(gt)à(0)e~ 



+ 



g\ cos(gt) + g\ 



+ 



9i92[cos(gt) - 1] 



& 2 (0)é 



-ÏLüt 



b2[t) = sM gt)à(0)e-^ + !] L ™„-*-t 



+ 
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92 cos(gí) + g\ 



6 2 (0)e 



-ïüjt 



(3a) 
(3b) 
(3c) 



with 9 = V9Ï+92- 

We assume that the cavity field is initially prepared in 
a superposition of two distinct coherent states \a>i) and 
\a 2 ) in the normalized form of 



toK0)) c = <7|ai>+£>|a 2 >. 



(4) 



|tf(0)> - (C|ai)+D|a 2 ))®|0)i|0) 2 . 



(5) 



The time dependent wave function |*(í)) of the whole 
system can usc the time evolution operator U(t) — 
e -i(Ht/h) and m j t i a i state |$(0)) to express as |*(í)) = 

U{t)\V(0)), that is 



|tf (i)> = U(t)(C\ ai ) + D\a 2 )) ® |0)i|0) 2 . 



(6) 



According to the definition of the coherent state, any co- 
herent state of the cavity field has the following form 



a) = exp^t (0) - a*a(0)]|0) 



(7) 



where (*) denotes complex conjugation. Then (6) can be 
rewritten as 

|*(í)) = U(t)[Cexp( ai à\0) - aífi(O)) 

+ Dexp(a 2 a t (0) - c^à(0))]|0)|0)i|0) 2 . (8) 

We may interpolate unit operator W{t)U{t) into Eq.(8) 
and consider the properties of the time evolution oper- 
ator W{t)OU(t) = 0(t) and U(t)\0) = |0), then Eq.(8) 
becomes as 



|tf(t)> = C\a lU (t)} <g> |ai«i(í))i <8> \a 1 v 2 {t)) 2 
+ D\a 2 u(t)) <g> |a 2 «i(í))i ® \a 2 v 2 (t)) 2 



with 



u(í) = cos(gt)e- iut , 

vUt) = -i^- sm(gt)e~ luJt , 
9 
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V2(t) 



-i— sin(gí)e 



(9) 

(10a) 
(10b) 

(10c) 



For convenience, we will drop the time dependence t from 
u(t), vi(t), and v 2 (t) in the following expressions. It 
is very clear that three subsystems become entangled 
with the time evolution except at some critical times 
t = niï/2g for integer n. But what we are interested 
in is the entanglcmcnt between two subsystems of exci- 
tons, so we trace over mode of the cavity field in the 
density matrix of the system given as \^ (t)} (t)\ to ob- 
tain the reduced density matrix p(t) = Tr(|*(í))(*(í)|) 
of excitons for any arbitrary time t. This partial trace 
operation enables us to describe the dynamics of the ex- 
citons in the coupled quantum dot system without any 
reference to the cavity field. Then the reduced density 
matrix p{t) is found as 

p(t) = [C*lW(|a2Ui)<aiUi|)i <g> {\a 2 v 2 )(a lV2 \) 2 + h.c] 
+ |C| 2 (|aii;i)(a!ií;i|)i <g> {\aiv 2 ) (aiv 2 \)2 
+ | J D| 2 (|a 2 wi)(a 2 wi|)i <g> (\a 2 v 2 )(a 2 v 2 \) 2 . (11) 
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Hcrc h.c. denotes complex conjugate. M 
product of two coherent states \a 2 u) and 



is an mncr 
ot\u), and 
. It is evi- 



M = exp [- cos 2 (5í)(i|ai| 2 + \\a 2 \ 2 - a*a 2 ) 
dent that unitary evolution of the wave function cannot 
change its purity, however tracing out the cavity field 
reduces the other subsystems into a statistical mixture. 



III. ENTANGLEMENT OF EXCITONIC STATES 

Concurrence is a useful tool to measure entanglemcnt 
bctween two subsystems of 2 x 2 biparte mixed and pure 
states. According to the references [16,17], we assume 
that 1 and 2 is a pair of qubits, and the density matrix 
of the pair is p\ 2 . Then the concurrence of the density 
matrix pi 2 is defined as: 



C12 = max{Xi — X 2 — A3 - A4, 0} 



(12) 



where Ai, A2, A3, and A4 given in decreasing order are 
the square roots of eigenvalues for matrices 



M12 = pi 2 (criy <8> o 2y )p\ 2 {o\ y ® <J 2y ) 
with Pauli matrix 



<Jly 



-i 

1 



<J 2 y 



-i 

1 



(13) 



(14) 



and (*) denoting complex conjugation in the Standard ba- 
sis {|00), |01), |10), |11)}, and a ly and a 2y are expressed 
in the same basis. The entanglement of formation is a 
monotonically increasing function of Ci 2 . C\ 2 = corre- 
sponds to an unentangled state, and C í2 = 1 corresponds 
to a maximally entangled state. 

First, we make a suitable transformation on Eq.(ll) 
so that we can define the qubit and discuss the en- 
tanglement between two modes of excitons using the 
concurrence. We can choose two different coherent 
states, which are nonorthogonal and form a super- 
complete set, as a basis to span a two-dimensional 



Hilbert space. But it is not suitable to define qubits 
using two-dimensional supcr-complcte set of basis. How- 
ever, one can always rebuild two orthogonal and nor- 
malized states as basis of the two-dimensional Hilbert 
space using original two coherent states [19]. In our 
system, we can choose two timc-dependent excitonic 
coherent states to define qubits. We define the zero 
qubit state as |0) m = \aiv m (t)}; and one qubit state as 
|l)m = {\a 2 v m {t )) - (aiv m {t)\a 2 v m (t))\ai Vm(t )))/N m (t) 

with N m (t) = ^l-\(a lVm (t)\a 2 v m (tW = x/Ï^PjJt) 
for exciton systems of mode one m = 1 and two m = 2. 
The above definition for qubits is equivalent to the def- 
inition in reference [20]. In this new basis, the reduced 
density matrix (11) can be rewritten into 

p = {cWKA^liKODx + /MíXloxoDx] ® 

[(JV 2 (í)|l)(0|) 2 +P 2 (í)(|0)(0|) 2 ]+/ l .c.} 
|D| 2 [(iVi(í)|l) + Pi(í)|0))(7V*(í)(l| + Pí(í)<0|)]i ® 
[(N 2 (t)\l) + P 2 (í)|0))(./V 2 *(í)(l| +P 2 (t)(0\)} 2 

+iq 2 (|o)(o|) 1 ®(|o)(o|)2. (15) 

For simplicity, we assume that two quantum dots are 
completely the same so that the interaction between cav- 
ity field and two quantum dots are equal to each other, 
that is gi = g 2 . We consider the cases where the cav- 
ity field is initially prepared in the even or odd coherent 
states. It has been shown that such mesoscopic coherent 
superpositions can be generated by experimentalists us- 
ing cavity QED [21] or a trapped ion [22]. For the cavity 
field initially prepared in the odd coherent state, Eq.(4) 
can be written as 



\odd) =N-(\a) - \-a)) 



(16) 



with normalization constant iV_ = (2 — 2e 2 H 2 ) 2. Un- 
der the above conditions, we can write N\ (t) = N 2 (t) = 
M = y/1 - P 2 {t) and P = Pi{t) = P 2 {t) = e -sín 2 (gt)\ a \\ 
Then the reduced density operator (15) can be rewritten 
in the matrix representation as follows 



R(p) = 



2(1 - e - 2 M 2 ) 



/ 1-2AÍP 2 + P A 
-MMP + MP 3 
-MMP + MP 3 
-MM 2 + M 2 P 2 



-MMP + MP 3 
M 2 P 2 
M 2 P 2 
M 3 P 



-MMP + MP 3 

M 2p 2 
M 2p 2 

M 3 P 



-MM 2 + M 2 P 2 
M 3 P 
M 3 P 
M 4 



(17) 



in the Standard basis {|00), |01), |10), |11)} where every 
basis vector such as 1 00) denotes state |0)i (g) |0) 2 - We 
can obtain the square roots of eigenvalues of the matrix 
(13) corresponding to the matrix (17) as 

(1 - e- 2sin2 (9*)l Q | 2 )(l + e -2cos 2 ( 3 í)H 2 ) 

^1 = — _oi„,i2x (18a) 



A 2 = 



2(l- e -2|a| 2 ) 

(1 - e — 2 sin 2 Cíï*)l"| 2 ) (1 _ e -2co S 2 ( ff t)|a| 2 -) 



2 (l_ e -2|a| 2 ) 



A 3 = A 4 = 0. 



(18b) 
(18c) 



Then the concurrence C of the prepared state with an 
initial odd coherent state of the cavity field is found as 

-2cos 2 ( ff t)|a| 2 ( '-| _ -2sin 2 ( S í)H 2 \ 

c- (é^n l - < w > 

We can also obtain the average photon number n of 
the cavity field as follows 
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Eqs. (19-20) show that two subsystems of the excitons 
becomes maximally cntanglcd for any time satisfying the 
condition gt = (n + for any intensity \a\ 2 of the 
coherent cavity field, and in this moment no photon of 
the cavity field can be detected. Fig. 1 shows the time 
evolution of (19) for two different intensity vàlues \a\ 2 of 
the cavity field. It is clearly seen that maximally entan- 
gled coherent states of excitons can be prepared by the 
cavity field of odd coherent states. The entanglement 
pcriodically reachcs its maximum. It is also found that 
when the excitonic states are in a maximally entangled 
state, the average photon number of the cavity field is 
zero (see Fig. 3). This enables us to monitor the status of 
the system during preparation of the entangled excitonic 
coherent states by detecting the photon number of the 
cavity field. 

If the cavity field is initially prepared in an even co- 
herent state 



jeuen) = N+(\a) 



a» 



(21) 



where the normalization constant is N + = (2 + 
2 , using the same procedure the concurrence 
C e is found as 



-2cos 2 (3t)|a| 2 ( '-| _ -2sin 2 ( 9 í)H 2 ï 

C e = „ U „ e ,, 9 , >-, (22) 



(1+6-2H 2 ) 



and the average number n of the cavity field is 



From Eq.(22), it is found that we can not obtain max- 
imally entangled states when the cavity field is initially 
in the even coherent state. But the analytical expres- 
sion (22) shows that when the coherent intensity \a\ of 
the cavity field tends to infinity, we can approximately 
obtain maximally entangled excitonic coherent state at 
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FIG. 1. Time evolution of the concurrence C is plotted as 
a function of — for (a) \a\ = 1 (dashed), (b) \a\ = 5 (sòlid). 
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FIG. 2. Time evolution of the concurrence C e is plotted as 
a function of gt/ir for (a) \a\ = 1 (dashed), (b) \a\ = 5 (sòlid). 
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FIG. 3. Time evolutions of the concurrence (sòlid curve) 
and average photon number n (dashed curve) are plotted as a 
function of gt/n for cavity field initially prepared with |a| — 1 
as odd coherent state (left) and even coherent state (right). 

time gt = (n + ^)n. In fact numerical results show that 
when the coherent intensity \a\ > 2, the concurrence 
approximately reaches its maximum at the same time 
gt = (n + |)7r. Fig. 2 shows the time evolution of the 
concurrence C e for two different \a\ vàlues of the cavity 
field. It is coneluded that a quasi-maximally entangled 
excitonic coherent states can be prepared using the initial 
even coherent state of the cavity field when the coherent 
intensity \a\ becomes intensive. 

In Fig. 3, a comparison of the degree of entanglement 
and its relation with the average number of photons (n) 
in the cavity for cavity fields initially prepared in odd 
and even coherent states is presented. It is seen that 
(i) Both the concurrence and n shows a periòdic evolu- 
tion, (ii) Concurrence reaches to its maximum value when 
n = for both cases, and (iii) for the cavity field initially 
prepared in odd coherent state, it is possible to obtain 
maximally entangled states for some evolution times, but 
this is not true for cavity field initially prepared in even 
coherent state. 



IV. EFFECT OF ENVIRONMENT ON THE 
ENTANGLEMENT OF EXCITONS 

In this section, we will focus our attention on how the 
environment with zero temperature affects the entangle- 
ment of the excitons. We assume that the energy of 
the system is dissipated only by the interaction of the 
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cxcitons in the quantum dots with the multi-modes of 
electromagnètic radiation. Our discussion is limited to 
these two completely equivalent quantum dots, that is, 
they have the same coupling constants with the cavity 
field, and also have the same dissipative dynamics. Undcr 
above considerations, we have the following Hamiltonian 
with the rotating-wave approximation as 



va— 1 k 

2 

+ Y K m.k(a m .kbln + a: m k b m ) 



m=l m,k 



where Hq is detcrmined by (1), and 



m.k( a "i,k 



dcnotc 



the multi-modes electromagnètic field operators with the 
freqüències u) m .k, an d m = 1(2) means the quantum dot 
one (two). Using the assumption g\ = g 2 = k — g/V%, 
and Ki y k — K 2,k, we can write the equations of motion for 
all operators as follows 



dà 
dt 

dk 



dt 

da m ,k 
dt 



= —iu) à — inb\ — inb 2 , 
- = -iuiob m - ina - i ^ K m ,ka m ,k, 



~ÍU) m ^k(^m^k i^m^kbm· 



(25a) 
(25b) 

(25c) 



We know that only the solution of the cavity field op- 
erator is enough to solve our problem. We apply the 
Laplacc transformation and the Wigner-Wcisskopf ap- 
proximation [23] to the above set of Eqs.(25a-25c). So 
the time-dependent solution of the cavity field operator 
with the zero-temperature environment can be obtained 
as 

a(t) = u'(t)a(0) + «í(t)6i(0) + v' 2 (t)b 2 (0) 

+ 5>i,fc(í)ai, fe (0) + 5> 2 ,fc(íHfc(0) ( 26 ) 



where 



«'(*) = 



e * 



í [ 7 (ei Aí -e-4 



)+íA(ei At + e-3 At )]e- í " t 



v[{t) =v' 2 (t) 



Í2A 

ÏA 



with A = y / 32n 2 — 7 2 , and 7 = 27re(u;o)|K(wo)| 2 where 
e(ui) is a distribution function of the multimode electro- 
magnètic radiation field. A small Lamb frequeney shift 
has been neglected in above Eqs.(27a-27b). We take 
the same steps as in the Sect. II, and obtain the time- 
dependent wave function of the whole system with cavity 
field initially in the odd mesoscopic coherent superposi- 
tion of states (16) and other subsystems in the vacuum 
states 



|*(t)) - C\au'(t)) ® I] \av' m (t)) ® ]J II I^WW) 

rn — 1 m— 1 k 

2 2 

-c\ - au'( t )) ® n i - av '^) ® n n i - ^w»' 



m—l 



m—1 k 



After tracing over the environment and the cavity field, 
we obtain the reduced density matrix of two subsystems 
of the cxcitons as 



( 24 ) P ' (t) - [|qW(|a 2 «í)(a>i|)i ® (\aW 2 ){ ai v' 2 \) 2 + h.c] 
+ |C| 2 (|an;í)<aiví|)i ® (\a lV ' 2 ) (a lV ' 2 \) 2 
+ |C| 2 (|a 2 t;i)<a2«i|)i ® (\a 2 v' 2 )(a 2 v' 2 \) 2 (29) 



with J\í' = e-Wfií-^vW 2 ) ^ where the conc jition 

E fe l«i, fc (t)l 2 + E fc l«2, fc (í)| 2 = i - KWI 2 r KWI 2 - 

|w2(í)| 2 , which comes from the commutation relation 
[a{t),a){t)] = 1, is used. Then we choose the basis as 
the section II, and use the same step to obtain the con- 
currence corresponding to the reduced density matrix op- 
erator (29) 



cí 



-2| Q | 2 (l-2|^| 2 )( 1 _ e -4|a| 2 |«' 1 | 2 - ) 



(1 



*l 2 ) 



(30) 



where v[ is determined by (27b). Whcn 7 = which 
means that the whole system has no any dissipation of 
energy, then we find that (30) returns to ideal case (19). 
Figs. 4 and 5 show that the cntanglcmcnt of the two 
coherent excitonic states is gradually reduced by the en- 
vironment with time evolution. If the decay rate 7 of 
the excitons is fixed, then the higher the intensity of the 
cavity field is, the faster the decay of entanglement is. In 
the same way, if the intensity of the cavity field is fixed, 
then for a system with larger decay rate 7 of excitons, 
the amount of entanglement reduces faster. 



(27a) 
(27b) 
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FIG. 4. Time evolution of the concurrence C is plotted as 
a function of gt/ir for -y/g — 0.01 when \a\ = 5 (sòlid), and 
\a\=2 (dashed). 
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FIG. 5. Time evolution of the concurrence C is plotted as 
a function of gt/-K for y/g = 0.01 (dashed), and y/g = 0.04 
(sòlid) when \a\ = 2. 



V. CONCLUSION 

The entanglcmcnt of formation for coherent states of 
the excitons in the system of two coupled quantum dots 
by cavity field is investigated. We find that a maximally 
entangled coherent state of two-modes excitons can be 
prepared by an odd coherent state of the cavity field. 
The status of this entangled excitonic coherent state can 
be probed by the measurement of the cavity field. No 
photon detection in the cavity implies the preparation 
of maximally entangled excitonic coherent states. In the 
case of photon detection of any number, it can be said 
that the excitons are not maximally entangled. 

The effects of the zero-tcmpcrature environment on the 
entanglement of excitons is also studied by using con- 
currence. It is found that the amount of entanglement 
between two excitonic coherent states is reduced by the 
environment with time evolution due to the loss of pho- 
tons of the cavity field to the environment. If leakage 
from the cavity to the environment is monitored by pho- 
ton counting detectors, an absenec of photo-detections 
implies that maximally entangled state of the excitons 
can be found in the coupled quantum dot system with 
time evolution. In case of any photon detection, we are 
certain that the system can not reach to a maximally 
entangled state. If the decay rate y of the excitons is 
given, for the system with higher cavity field intensity, 
the decrease in the amount of entanglement is faster. It 
is also understood that if the coherent intensity of the 
cavity field is given, then the large decay rate 7 of exci- 
tons corresponds to a fast reduction in entanglement of 
the system. 

To have an idea on the characteristic time-scales of 
the dynamics of entanglement in our model, we can take 
hg ~ 0.5 meV which may be possible to obtain in a sys- 
tem of large area quantum dot with cavity [7]. In that 
case, for an ideal cavity with initial field prepared as an 
odd coherent state, the first maximally entangled state 
could be observed after a time evolution of 2 ps and this 
would repeat itself periodically with a period of 4 ps. For 
a practical cavity, the finite cavity lifetime may be taken 



as 7 _1 ~ 10 ps [7]. For these g and 7 vàlues, we can find 
■y/g ~ 0.13. With \a\ = 2 (|a| = 1), the maximum value 
for concurrence would be ~ 0.46 (<~ 0.8)which oceurs at 
t ~ 2 ps. After a time evolution of 4 ps, the second peak 
for the concurrence would appear with a value of ~ 0.17 
(<~ 0.53). It is clearly seen that for a practical scheme 
the main parameters which would affect the dynamics of 
the system are the fast rate of the photon loss from the 
cavity and intensity of the initial cavity field. 
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